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Abstract. We report on diffuse X-ray emission associated with the nearby solitary millisecond pulsar PSR J2 124-3358 detected 
with XMM-Newton and Chandra. The emission extends from the pulsar to the northwest by ~ 0.5 arcmin. The spectrum of 
the nebular emission can be modeled with a power-law of photon index 2.2 + 0.4, in line with the emission originating from 
accelerated particles in the post shock flow. For PSR J0437-4715, PSR J0030+0451 and PSR J1024-0719, which all have 
spin parameters comparable to that of PSR J2124-3358, no diffuse emission is detected down to a 3 - a limiting flux of 
~ 4 - 7 x 10~ 15 ergs s _1 cnr 2 . 

Key words, pulsars: individual (PSR J2124-3358, J0437-4715, J0030+0451, J 1024-07 19)— -stars: neutron— X-rays: stars 



1. Introduction 

Diffuse plerionic emission is often a special marker of young 
and powerful pulsars if their spin-down energy is in excess 
of ~ 10 36 erg/s. Extended diffuse emission associated with 
older and less powerful pulsars, however, is rare and so far 
only seen in Geminga (Caraveo et al. 2003), PSR B1929+10 
(Becker et al. 2005), PSR B2224+65 (Romani et al. 1997; 
Chatterjee & Cordes 2002) and the millisecond pulsars (MSPs) 
PSR B 1957+20 (Stappers et al. 2003), J0437-4715 (BeU et 
al. 1995), and J2124-3358 (Gaensler, Jones & Stappers 2002). 
While all three MSPs have bright bow-shock nebulae detected 
in H a , diffuse X-ray emission associated with it could be de- 
tected only from the black widow pulsar PSR B 1957 +20 in a 
recent Chandra observation (Stappers et al. 2003). Comparable 
deep observations of almost all X-ray bright MSPs, however, 
have been performed by XMM-Newton in previous years. 
These XMM-Newton data thus provides us a valuable data base 
to search for diffuse and extended emission components from 
these MSPs as well. 

In this letter we report on a search for diffuse X-ray emis- 
sion to be associated with the MSPs PSR J2124-3358, J0437- 
4715, J0030+0451 and J1024-0719. All these pulsars have 
comparable spin parameters (cf. Table 1) so that differences in 
their emission properties are most likely cause by differences 
in the pulsar-ISM interaction/local environment rather than by 
differences in their total energy output. 

PSR J2124-3358 and J1024-0719 were both discovered 
during the Parkes 436 MHz survey of the southern sky (Bailes 
et al. 1997). PSR J0030+0451 was discovered at 430 MHz 
during the Arecibo Drift Scan Search (Somer 2000) and in- 
dependently in the Bologna sub-millisecond pulsar survey 
(D'Amico 2000), whereas PSR J0437-4715 was discovered in 



the Parkes southern sky survey (Johnston et al. 1993). PSR 
J0437-4715 was the first MSP of which pulsed X-ray emission 
was detected (Becker & Triimper 1993). X-ray emission from 
PSR J2124-3358 and J1024-0719 was reported by Becker & 
Triimper (1998, 1999) in ROSAT HRI data whereas the X- 
ray counterpart of PSR J0030+0451 was discovered in the 
final ROSAT PSPC observation (Becker et al. 2000). A H a 
bow shock is seen around PSR J0437-4715 (Bell et al. 1995). 
However, its X-ray counterpart was not detected by ROSAT 
and Chandra (Becker & Triimper 1999, Zavlin et al. 2002). 
Gaensler, Jones & Stappers (2002) discovered an //^-emitting 
bow shock nebula around PSR J2124-3358. This bow shock 
is very broad and highly asymmetric about the direction of 
the pulsar's proper motion. The asymmetric shape might be 
caused by a significant density gradient in the ISM, bulk flow 
of ambient gas and/or anisotropics in the pulsar's relativistic 
wind (Gaensler, Jones & Stappers 2002). Observation of PSR 
J0030+0451, performed with the ESO NTT in La Sila, did 
not detect diffuse H a emission associated with it (A. Pellizzoni 
priv. com.). 



2. Observations and data analysis 

The pulsars PSR J2124-3358, J0030+0451, J1024-0719 and 
J0437-4715 were observed by XMM-Newton in 2002 April 
14-15, in 2001 June 19-20, in 2003 December 2 and in 2002 
October 9, respectively. In all observations the MOS1/2 cam- 
eras were operated in full-frame mode while the PN camera 
was setup to work in the fast-timing mode. The PN fast-timing 
mode provides only limited spatial information so that these 
data are not used in the present analysis. To block stray light 
and optical leakage from bright foreground stars we operated 
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Table 1. Pulsar Parameters of PSRs J0030+0451, J2124-3358, J1024-0719 and J0437-4715 (from Manchester et al. 2005) 



Pulsar 


PSR J0030+0451 


PSR J2124-3358 


PSR J1024-0719 


PSR J0437-4715 


Right Ascension (J2000) 


00 h 30 m 27.432 s 


21 h 24 m 43.862 s 


10 h 24 m 38.700 s 


04 h 37 m 15.787 s 


Declination (J2000) 


+04° 51' 39.7" 


-33° 58' 44.257" 


-07° 19' 18.915" 


-47° 15' 08.462" 


Pulsar Period, P (ms) 


4.865453207369 


4.9311148591481 


5.1622045539088 


5.7574518310720 


Period derivative P (1CT 20 s s~') 


1.0 


1.33 


1.85 


1.87 


Age (10 9 yrs) 


7.71 


5.86 


4.41 


4.89 


Surface dipole magnetic field (10 s G) 


2.23 


2.60 


3.13 


3.32 


Epoch of Period (MJD) 


50984.4 


50288.0 


51018.0 


51194.0 


Dispersion Measure (pc cm 4 ) 


4.3328 


4.6152 


6.491 


2.6469 


Dispersion based distance (pc) 


230 


250 


350 


140 


Spin-down Luminosity (10 33 ) ergs s~ 2 


3.43 


4.38 


5.3 


3.87 



the MOS1/2 cameras with the thin filters during the observa- 
tions of PSRs J0030+0451, J 1024-07 19 and J0437-4715 while 
PSR J2 124-3358 was observed by using a medium filter. Events 
were selected for the energy range 0.25 - 10 keV and stan- 
dard correction procedures were applied to reduce the data 
(e.g. Becker et al. 2005). Filtering the data for times of exces- 
sive background from soft proton flares was done using con- 
servative thresholds in view of the planned search for diffuse 
and extended X-ray emission. The effective exposure time af- 
ter data reduction turns out to be 40 ksec for PSR J2124-3358, 
14 ksec for J0030+0451, 65 ksec for J1024-0719 and 55 ksec 
forJ0437-4715. 

A vignetting corrected image of the field of PSR J2124- 
3358 as seen by the XMM-Newton's MOS1/2 CCDs is shown 
in Figure la. The binning factor in this image is 6 arcsec. 
Adaptive smoothing with a Gaussian kernel of cr < 4 pixels 
has been applied to the image. X-ray contours are calculated 
and overlaid on the image. The contour lines are at the levels 
of (4.2,5.2,7.6, 13,28,63) x 10~ 6 cts s" 1 arcsec- 2 . It can be 
seen that the X-ray source which is coincident with the pul- 
sar position has an asymmetric source structure of ~ 0.5 ar- 
cmin extent, with its orientation to the northwest. Systematic 
effects which could cause an adequate distortion of the instru- 
ment's point spread function (PSF) are not known for XMM- 
Newton. We are therefore prompted to interpret this elongated 
structure in terms of a pulsar X-ray trail. The net count rate 
of the diffuse emission in a 30 x 35 arcsec box centered at 
RA(J2000)=21 h 24 m 42.69 s , Dec=-33°58'16.06" and oriented 
along the extended feature is estimated to be (1.25±0.16)xl0~ 3 
cts s _1 . For this estimate we determined the background from 
a low count region close to PSR J2124-3358. For compari- 
son, the pulsar emission is estimated to have a net count rate 
of (1.07 + 0.01) x 10 -2 cts s~ : in a circle of 18 arcsec radius 
centered at the pulsar position. The signal-to-noise of this elon- 
gated feature is ~ 4 in the energy range 0.25 - 5 keV. 

PSR J2124-3358 was observed by Chandra with the back- 
illuminated ACIS-S3 chip in the focus of the XRT The ob- 
servation took place on 2004 December 19-20 for an expo- 
sure time of ~ 30 ksec. An image made from this data with 
0.5 arcsec binning and adaptive smoothing applied (using a 
Gaussian kernel with cr < 1.5 pixel) is shown in Figure lb. 
Arc-like diffuse emission which is within the pulsar's H a neb- 
ula is clearly detected. Its net count rate is estimated to be 
( 1 .28 + 0.23) x 1 0~ 3 cts s" 1 in a box of 1 x 20 arcsec centered at 
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Fig.l. a) XMM-Newton MOS1/2 image of PSR J2124-3358 
with overlaid contours. The pulsar proper motion is indi- 
cated by an arrow. The position of bright stars located in the 
1.5 arcmin neighborhood of PSR J2124-3358 are indicated, 
b) PSR J2124-3358 as seen by the ACIS detector aboard 
Chandra. 
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Table 2. Identifications of the stars around the X-ray trail of 
PSRJ2124-3358. 



Stars " 


USNO catalogue ID 


Magnitude: R 


Magnitude: B 


A 


U052543254607 


17.1 


18.7 


B 


U052543254423 


18.0 


20.2 


C 


U052543253997 


16.7 


17.4 


D 


U052543253910 


13.5 


15.5 



a. See Figure la. 



RA(J2000)=21 h 24 m 43.17 s , Dec=-33°58'33.53". With the aid 
of PIMMS, the pileup fraction of the ACIS-S3 data was found 
to be less than 0.2%, i.e. is negligible. The signal-to-noise of 
this feature is ~ 5 in the energy range of 0.3 - 8 keV. 

Since the PSF width of XMM-Newton is about 10 times 
that of Chandra, it blurred most of the detailed structure seen in 
the Chandra data. However, it should be noted that the overall 
direction of the feature in the Chandra image is consistent with 
the orientation of the trail detected by XMM-Newton. 

In order to investigate a possible contribution to the dif- 
fuse X-ray emission by nearby stars we investigated the 
Digitized Sky Survey data (DSS) for the sky region around 
PSR J2124-3358. There are four field stars (labeled as A, B, C 
and D in Figure 1 a and Table 2) in the 1 .5 arcmin neighborhood 
of PSR J2124-3358. From the USNO catalogue, we identified 
the magnitudes of these stars. None of them is found to match 
the position of the diffuse elongated X-ray structure seen in the 
XMM-Newton and Chandra data. It can be seen in Figure la 
that the positions of stars B and C coincide with two faint X- 
ray sources which are disconnected with the trail emission of 
PSRJ2124-3358 though. 

To examine the nature of the faint X-ray nebular emission 
by a spectral analysis we make use of the XMM-Newton data 
and extract the energy spectrum from within a 30 X 35 arcsec 
box. Using a box rather than a circular selection region allows 
to avoid the emission from the pulsar and excludes any poten- 
tial contamination from the field stars B and C. However, we 
estimate that the wing of the XMM-Newton PSF centered at the 
pulsar position still contributes ~ 20% to the total counts inside 
the box. The background spectrum was extracted from a source 
free region near to PSR J2124-3358. In total, 92 and 67 source 
counts are available from the trail in the MOS1 and MOS2 
detectors, respectively. Response files were computed by us- 
ing the XMMSAS tasks RMFGEN and ARFGEN. The spectra 
were dynamically binned so as to have at least 10 counts per 
bin. 

In the Chandra data we selected the energy spectrum of the 
diffuse emission from a box of size 10x20 arcsec. Owing to the 
narrow PSF of Chandra the contamination of pulsar emission in 
this box is negligible. The background spectrum was extracted 
from a low count region near to the diffuse feature. In total 46 
source counts are contributed from the Chandra data. Response 
files were computed by using tools MKRMF and MKARF of 
CIAO. The spectrum was binned to have at least 8 counts per 
bin. The degradation of quantum efficiency of ACIS was cor- 
rected by applying the XSPEC model ACISABS. 

We hypothesize that the diffused emission should be orig- 
inated from the interaction of pulsar wind and the ISM. 




0.5 1 2 

channel energy (keV) 



Fig. 2. Energy spectrum of the X-ray trail of PSR J2124-3358 
as observed with MOS1/2 and ACIS-S3 detectors and simulta- 
neously fitted to an absorbed power-law model (upper panel) 
and contribution to the^- 2 fit statistic (lower panel). 

Synchrotron radiation from the ultra-relativistic electrons is 
generally believed to be the emission mechanism of the pul- 
sar wind nebula, which then is characterized by a power-law 
spectrum. With a view to test this hypothesis, we fitted an ab- 
sorbed power-law model to the nebular spectrum with XSPEC 
11.3.1 in the 0.25 - 10 keV energy range. With a column den- 
sity of 5 x 10 20 cirT 2 as obtained from spectral fits to the pul- 
sar emission, we found that the model describes the observed 
spectrum fairly well (xl = 0.79 for 26 D.O.F.). The photon 
index is a — 2.2 + 0.4 and the normalization at 1 keV is 
(2.94 + 0.48) x 10~ 6 photons keV" 1 cirT 2 s _1 (1 - cr error for 
1 parameter in interest). In view of the low photon statistic we 
tested for a possible dependence of the fitted model parameters 
against the background spectrum. All deviations found in re- 
peating the fits with different background spectra were within 
the 1 — cr confidence interval quoted above. The unabsorbed 
fluxes and luminosities deduced for the best fit model param- 
eters and the energy ranges 0.1-2.4 keV and 0.5-10 keV are 
fx = 1.8 X 10~ 14 ergs s _1 cirT 2 , L x - 13 x 10 29 ergs s _1 and 
fx= 1.2 x 10~ 14 ergs s" 1 cirT 2 , L x = 8.9 x 10 28 ergs s _1 , re- 
spectively. The best fitting spectral model is displayed in Figure 
2. 

Motivated by the extended source structure seen in 
PSR J2124-3358 we have also inspected the XMM-Newton 
data of PSRs J0030+045 1 , J 1 024-07 1 9 and J0437-47 1 5 in or- 
der to search for extended tail-like X-ray emission. There is 
no evidence for any extended emission around these pulsars. 
In order to deduce an upper limit for an X-ray trail in these 
sources we have estimated the net count rates for the regions 
of 30 x 35 arcsec right behind the corresponding directions of 
each pulsar's proper motion. The net count rate for any dif- 
fuse emission in these regions is (3.1 + 1.5) x 10~ 4 cts s _1 , 
(3.1 ± 0.7) x 10~ 4 cts s- 1 and (2.5 + 0.9) x 10~ 4 cts s" 1 for 
PSRs J0030+0451, J 1024-07 19 and J0437-4715, respectively. 
The quoted errors are 1 - cr confidence intervals. Assuming a 
Crab-like spectrum (i.e. a = 2) and taking the values of column 
density inferred from the pulsar spectral fits (1 x 10 20 cm -2 for 
J0030+0451, 2 x 10 20 cirT 2 for J1024-0719 and 3 x 10 19 crrT 2 
for J0437-4715), we calculated the 3 - cr limiting fluxes from 
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the above countrates. For PSRs J0030+0451, J1024-0719 and 
J0437-4715, the limiting fluxes in 0.5-10 keV are estimated to 
be 6.9 x 10~ 15 ergs s -1 cirT 2 , 4.9 x 10~ 15 ergs s _1 cirT 2 and 
4.2 x 10~ 15 ergs s _1 cirT 2 , respectively. 

3. Discussion & Conclusion 

Basic characteristics of the X-ray bow shocks detected 
around rotation-powered pulsars are summarized in Table 3. 
PSR J2124-3358 is a new candidate to this group. Its tail is 
found to have an extent of ~ 0.5 arcmin. Adopting the distance 
of ~ 250 pc, the tail has a length of I ~ 1.1 x 10 17 cm. For 
the pulsar's proper motion velocity of 58 km s _1 (Manchester 
et al. 2005), the timescale, ffl ow , for the passage of the pulsar 
over the length of its X-ray trail is estimated to be ~ 600 yrs. 
According to the discussion in Becker et al. (2005) on the trail 
emission of PSR B1929+10 we estimate the magnetic field 
in the shocked region by assuming ffl ow to be comparable to 
the electron lifetime of the synchrotron emission. This yields 
~ 30yuG for the inferred magnetic field strength in the emitting 
region. The magnetic field strength in the ISM is estimated to 
be ~ 2 - 6fiG (cf. Beck et al. 2003 and references therein). 
Taking into account that the magnetic field in the termination 
shock might be compressed (e.g. Kennel & Coroniti 1984), our 
estimation is approximately consistent if the compression fac- 
tor is ~ 7. 

Following Becker et al. (2005), we apply a simple one 
zone model (Chevalier 2000) to estimate the spectral behav- 
ior and the X-ray luminosity of the nebular emission. The X- 
ray luminosity and spectral index depend on the inequality be- 
tween the characteristic observed frequency v° bs and the elec- 
tron synchrotron cooling frequency v c which is estimated to be 
1.6 x 10 17 Hz. Since in general v* s > v c , we concluded that 
the emission is in a fast cooling regime. Electrons with the en- 
ergy distribution, N(~y) oc y~ p , are able to radiate their energy 
in the trail with photon index a = (p + 2)/2. The index p due 
to shock acceleration typically lies between 2 and 3 (cf. Cheng, 
Taam, & Wang 2004 and references therein). Taking p = 2.35 
yields a th ^2.2 which is in accordance with the result from the 
observed value Q' obs = 2.2 + 0.4. Assuming the energy equipar- 
tition between the electron and proton (Cheng, Taam, & Wang 
2004), we take the fractional energy density of electron e e to 
be ~ 0.5 and the fractional energy density of the magnetic field 
€b to be ~ 0.01. Assuming a number density of ISM to be 1 
cm -3 , the distance of the shock from the pulsar is estimated to 
be ~ 3.6 x 10 16 cm. With these estimates, the calculated lumi- 
nosity, vL v , is given as ~ 10 29 ergs s _1 which is well consistent 
with the observed values of 1.3 x 10 29 ergs s _1 (0.1-2.4 keV) 
and 8.9 x 10 28 ergs s" 1 (0.5-10 keV). 

Although the general properties of the X-ray trail in 
PSR J2124-3358 are not in contradiction with properties ob- 
served in other pulsars there are still ambiguities which are not 
completely resolved. First, one should notice that the trail is 
misaligned with the direction of the pulsar's proper motion. As 
reported by Gaensler, Jones, & Stappers (2002), the head of the 
H a bow shock is found to be highly asymmetric about the pul- 
sar's velocity vector with the apparent nebular symmetry axis 
deviated from the velocity vector by ~ 30°. Even though the 



misalignment of the X-ray trail seems to agree with the asym- 
metry of the H a nebula, deeper observations by XMM-Newton 
and Chandra are required in order to constrain the physical 
properties of this interesting nebula in higher detail. 
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PSR 


Other Name/ 


E 


L%" (0.5-10keV) a 


T spin-down 


H„ 


X-ray 


Radio 


Reference 




associated SNR 


erg s -1 


erg s -1 


years 


bow shock 


trail 


nebula 




Young pulsars 


J0537-6910 


LMC: N157B 


4.8 x 10 38 


3.5 x 10 36 


5.0 x 10 3 




d 




1,2 


B 1757-24 


Duck 


2.6 x 10 36 


4.6 x 10 32 


1.6 x 10 4 




d 


d 


3 


B1853+01 


W44 


4.3 x 10 35 


5.7 x 10 32 


2.0 x 10 4 




d 


d 


4 


J1747-2958 


Mouse 


2.5 x 10 36 


5.3 x 10 34 


2.6 x 10 4 




d 


d 


5 


B1951+32 


CTB80 


3.7 x 10 36 


2.5 x 10 33 


1.1 x 10 5 


d 


d 


d 


6,7,8 


Middle age and old pulsars 


B0740-28 




1.4 x 10 35 




1.6 x 10 5 


d 






9 


B0633+17 


Geminga 


3.2 x 10 34 


8.9 x 10 28 


3.4 x 10 5 




d 




10 


B2224+65 


Guitar 


1.2 x 10 33 


3.7 x 10 31 


1.1 x 10 6 


d 


? 




11,12 


B 1929+ 10 




3.9 x 10 33 


8.3 x 10 29 


3.1 x 10 6 




d 


? 


13 


Millisecond pulsars 


B 1957+20 




1.1 x 10 35 


1.9 x 10 31 


2.2 x 10 9 


d 


d 




14 


J0437-4715 




3.9 x 10 33 




4.9 x 10 9 


d 






15,16 


J2124-3358 




4.4 x 10 33 


10 29 


5.9 x 10 9 


d 


d 




17,18 



a. The X-ray luminosities of the diffuse emission are taken from the corresponding reference and recalculate into energy band of 0.5-10 keV 
for easy comparison. References: (1) Wang & Gotthelf 1998; (2) Wang et al. 2001; (3) Kaspi et al. 2001; (4) Petre et al. 2002; (5) Gaensler 
et al. 2004; (6) Migliazzo et al. 2002; (7) Moon et al. 2004; (8) Li et al. 2005; (9) Jones et al. 2002; (10) Caraveo et al. 2003; (11) Romani et 
al. 1997; (12) Chatterjee & Cordes 2002; (13) Becker et al. 2005; (14) Stappers et al. 2003; (15) Bell et al. 1995; (16) Zavlin et al. 2002; (17) 
Gaensler et al. 2002; (18) This letter; and Manchester et al. 2005 otherwise. 



